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ABSTRACT 

In this paper, we propose a novel channel estimation tech- 
nique based on 2D spread pilots. The merits of this tech- 
nique are its simplicity, its flexibility regarding the transmis- 
sion scenarios, and the spectral efficiency gain obtained com- 
pared to the classical pilot based estimation schemes used in 
DVB standards. We derive the analytical expression of the 
mean square error of the estimator and show it is a function 
of the autocorrelation of the channel in both time and fre- 
quency domains. The performance evaluated over a realistic 
channel model shows the efficiency of this technique which 
turns out to be a promising channel estimation for the future 
mobile video broadcasting systems. 

1. INTRODUCTION 

Orthogonal frequency division multiplexing (OFDM) has been 
widely adopted in most of the digital video broadcasting stan- 
dards as DVB-T [1], DMB-T, ISDB-T. This success is due to 
its robustness to frequency selective fading and to the simplic- 
ity of the equalization function of the receiver. Indeed, by im- 
plementing inverse fast Fourier transform (IFFT) at the trans- 
mitter and FFT at the receiver, OFDM splits the single chan- 
nel into multiple, parallel intersymbol interference (ISI) free 
subchannels. Therefore, each subchannel, also called subcar- 
rier, can be easily equalized by only one coefficient. 

To equalize the signal, the receiver needs to estimate the 
channel frequency response for each subcarrier. In the DVB- 
T standard, some subcarriers are used as pilots and interpo- 
lating filtering techniques are applied to obtain the channel 
response for any subcarrier. Nevertheless, these pilots reduce 
the spectral efficiency of the system. To limit this problem, 
we propose to add a two dimensions (2D) linear precoding 
(LP) function before the OFDM modulation. The basic idea 
is to dedicate one of the precoding sequences to transmit a 
so-called spread pilot [2] that will be used for the channel es- 
timation. The merits of this channel estimation technique are 
not only due to the resource conservation possibility, but also 
to the flexibility offered by the adjustable time and frequency 
spreading lengths. In addition, note that the precoding com- 
ponent can be exploited to reduce the peak-to-average ratio 



(PAPR) of the OFDM system [3], or to perform the frequency 
synchronisation. 

The contribution of this article is twofold. First, a gen- 
eral framework is proposed to describe the 2D precoding tech- 
nique used for channel estimation. Secondly, exploiting some 
properties of random matrix and free probability theories [4] 
[5], an analytical study of the proposed estimation method is 
presented. 

The article is organized as follows. In section 2, we present 
the principles of 2D LP OFDM, and detail the channel estima- 
tion technique using the spread pilots. In section 3, we anal- 
yse the theoretical performance of this channel estimation by 
developing the analytical expression of its mean square error 
(MSB). Then, simulation results in terms of MSB and bit error 
rate (BBR) are presented and discussed in section 4. Conclud- 
ing remarks are given in section 5. 

2. SYSTEM DESCRIPTION 
2.1. 2D LP OFDM 

Fig. 1 exhibits the proposed 2D LP OFDM system exploit- 
ing the spread pilot channel estimation technique. First of all, 
data bits are encoded, interleaved and converted to complex 
symbols Xjn,s [«]■ These data symbols are assumed to have 
zero mean and unit variance. They are interleaved before be- 
ing precoded by a Walsh-Hadamard (WH) sequence C; of L 
chips, with < i < L — 2" and n e N. The chips obtained 
are mapped over a subset of L = Lf.Lf subcarriers, with Lt 
and L f the time and frequency spreading factors respectively. 

The first Lt chips are allocated in the time direction. The 
next blocks of Lt chips are allocated identically on adjacent 
subcarriers as illustrated in Fig. 2. Therefore, the 2D chip 
mapping follows a zigzag in time. Let us define a frame as a 
set of Lt adjacent OFDM symbols, and a sub-band as a set of 
Lf adjacent subcarriers. In order to distinguish the different 
subsets of subcarriers, we define m and s the indexes referring 
to the frame and the sub-band respectively, with < s < 
5* — 1. Given these notations, each chip y„i^s [n, q] represents 
the complex symbol transmitted on the nth subcarrier during 
the qth OFDM symbol of the subset of subcarriers [m, s], with 
0<n<Lf — 1 and < q < Lt — I. Hence, the transmitted 
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Fig. 1. 2D LP OFDM transmitter and receiver based on 
spread pilot channel estimation technique 
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Fig. 2. 2D chip mapping scheme 



signal on a subset of subcarriers [to, s] writes: 

where Xm.s = [xm,s [0] . . . a::™,^ [«]••• Xm,s [L - 1]]^ is the [L x 1] 
complex symbol vector, P = diag { . . . ^/Pi ■ ■ ■ V-^i-i} is 
a [L X L] diagonal matrix where Pi is the power assigned to 
symbol x„i,s [i], and C = [co . . . . . . Cl-i] is the WH pre- 
coding matrix whose zth column corresponds to ?th precoding 
sequence Ci = [0,0] . . .Cj [n,q] . ..Ci[Lf-l,Lt-l]]'^. We as- 
sume normalized precoding sequences, i.e. [n, q] = i"^- 
Since the 2D chip mapping applied follows a zigzag in time, 
Ci [n, q] is the {n x Lt + g)th chip of the ith precoding se- 
quence c,. 

2.2. Spread pilot channel estimation principles 

Inspired by pilot embedded techniques [6], channel estima- 
tion based on spread pilots consists of transmitting low level 
pilot-sequences concurrently with the data. In order to reduce 
the cross-interferences between pilots and data, the idea is to 
select a pilot sequence which is orthogonal with the data se- 
quences. This is obtained by allocating one of the WH orthog- 
onal sequences Cp to the pilots on every subset of subcarriers. 
Let Hm,s be the [L x L] diagonal matrix of the channel coeffi- 
cients associated to a given subset of subcarriers [to, s]. After 
OFDM demodulation and 2D chip de-mapping, the received 
signal can be expressed as: 

2^m,s — ^^m,s^ m,s H" ^m^s (^) 

where Wr„,s = [wm,s [0, 0] . . . Wm,s [n, q]... Wr»,s[-L/-1, -L*- 
1]]^ is the additive white Gaussian noise (AWGN) vector hav- 
ing zero mean and variance cr^ = [n, olf^- 

At the reception, the de-precoding function is processed 
before equalization. Therefore, an average channel coefficient 



Havg [m, s] is estimated by subset of subcarriers. It is ob- 
tained by de-precoding the signal received Z^.s by the pilot 
precoding sequence Cp and then dividing by the pilot symbol 

Xm]s = \/PpXm,s [p] loiown by the receiver: 

^avg [^; ^] ~ (p) '^'in,s 
Xm,s 

Xm.s 

Let US define = [co . . . Ci^p . . . Cl_i] the [L x (L — 1)] 
data precoding matrix, = diag { ^/Po . . . y^Pi^p ■ ■ ■ ^y Pl-i} 
the [(L — 1) X (L — 1)] diagonal matrix which entries are the 
powers assigned to the data symbols, and x^.s = [a;m,s[0] . . . 
Xm,3 [i=^p]... Xm,s [L — 1]]^ the [{L — 1) x 1] data symbols 
vector. Given these notations, (3) can be rewritten as: 

Havg [m, s] 

- ^— fc^H c r^f) -I- c^H C P x^") + c^w 1 

= Havg [m, s] + SI [m, s] + w' (4) 

The first term Havg [w, s] is the average channel response 
globally experienced by the subset of subcarriers [m, s]. The 
second term represents the self-interference (SI). It results 
from the loss of orthogonality between the precoding sequences 
caused by the variance of the channel coefficients over the 
subset of subcarriers. In the sequel, we propose to analyse its 
variance. 

3. THEORETICAL PERFORMANCE OF THE 
ESTIMATOR 

In order to analyse the theoretical performance of the pro- 
posed estimator, we evaluate its MSB under the assumption 



of a wide-sense stationary uncorrelated scattering (WSSUS) 
channel. 



MSE [m, s] = £■ I Havg [m, s] - Havg [m, s] 

= E[\Sl[m,s]\'}+E[\w'f} (5) 
First, let us compute the SI variance: 

E{\SI[m,s]\'} = y^E [c^H^,,CuKcX,s<^p} (6) 

where P^, = P„Pf = diag {Pa... P^^p . . . Pl-i]. Actually, 
(6) cannot be analyzed practically due to its complexity. Ap- 
plying some properties of random matrix and free probability 
theories [4] [5] which is stated in Appendix, a new SI variance 
formula can be derived: 

e[\S\ [m, s]\^} = ^E H„,« (/ - Cpcf ) H^,,Cp} 



ttE { 



(7) 

The expectation of A is the average power of the channel co- 
efficients on the subset of subcarriers [m, s]. Assuming that 
the channel coefficients are normahzed, its value is one: 

. , Lf-lLt-l 

E Utr {Hm,sllZ,s)\ = 2 E T.^{\^rn,s[n,q]\^} 

n=0 g=0 

= 1 (8) 

The expectation of i? is a function of the autocorrelation of 
the channel Rhh (An, Aq) whose expression is developed 
in Appendix. Indeed, it can be written: 



Lf-lLt-lLt-lLt-l 



E{tr(H^,.)tr(H^,,)} = ^ E E E^«»(^"'^9) 

n=0 q=0 n'=0 q' =0 

(9) 

where An = n — n' and Aq = q — q'. Note that the autocor- 
relation function of the channel does not depend on the subset 
of subcarriers since the channel is WSSUS. By combining (8) 
and (9), the SI variance expression (7) can be expressed as: 
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Now, let us compute the noise variance: 



(11) 



Finally, by combining the expressions of the SI variance 
(10) and ttie noise variance (11), the MSE (5) writes: 

/ Lf-l Lf-l Lf-l Lt-l ^ 

MSE=-1 1--!^ E E E ^««(^"'^«)+'^» 



Pi" 

\ n=0 q=0 n'=0 q'=0 



(12) 



(10) 



The analytical expression of the MSE of our estimator de- 
pends on the pilot power, the autocorrelation function of the 
channel and the noise variance. The autocorrelation of the 
channel (18) is a function of both the coherence bandwidth 
and the coherence time. We can then expect that the proposed 
estimator will be all the more efficient than the channel co- 
efficients will be highly correlated within each subset of sub- 
carriers. One can actually check that if the channel is flat over 
a subset of subcarriers, then the SI (10) is null. Therefore, 
it is important to optimize the time and frequency spreading 
lengths, Lt and Lf, according to the transmission scenario. 

4. SIMULATION RESULTS 

In this section, we analyse the performance of the proposed 
2D LP OFDM system compared to the DVB-T standard un- 
der the COST207 Typical Urban 6 paths (TU6) channel model 
depicted in Table 1 with different mobile speeds. We de- 
fine the parameter as the product between the maximum 
Doppler frequency Jd and the total OFDM symbol duration 
ToFDM- Table 2 gives the simulation parameters and the useful 
bit rates of the DVB-T system and the proposed system. 

In the proposed system, only one spread pilot symbol is 
used over L > 16, whereas the DVB-T system uses one pilot 
subcarrier over twelve. Therefore, a gain in terms of spec- 
tral efficiency and useful bit rates are obtained compared to 
the DVB-T system. These gains are all the higher than the 
spreading factor L is high. Nevertheless, an increase of the 
spreading length produces a higher SI value. Consequently, a 
trade-off has to be made between the gain in term of spectral 
efficiency and the performance of the channel estimation. 

Fig. 3 depicts the estimator performance in term of MSE 
for QPSK data symbols, different mobile speeds and differ- 
ent spreading factors. The curves represent the MSE obtained 
with the analytical expression (12), and the markers those ob- 
tained by simulation. We note that the MSE measured by 
simulation are really closed to those predicted with the MSE 
formula. This validates the analytical development made in 
section 3. We note that beyond a given ratio of the energy per 
bit to the noise spectral density (^), the MSE reaches a floor 
which is easily interpreted as being due to the SI (5). 



Table 1. Profile of TU6 channel 
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Table 2. Simulation Parameters and Useful Bit Rates 



R fi n H w 1 H t Vi 


8 MHz 




2048 samples 


Guard Interval size 


512 samples (64 jis) 


iihIiA/T c\7Tnnr»l Hiirnfinn (' 1 i-ir-r-,* » 1 
WFl-'iVl syiiiUL'l ULIlaLIUli l^J-OFDM) 




Rate of convolutional code 


1/2 using (133, 171)^ 


Constellations 


QPSK and 16QAM 


Carrier frequency 


500 MHz 


Mobile Speeds 


20 km/hand 120km/h 


Maximum Doppler frequencies {fo) 


9.3 Hz and 55.6 Hz 


(3 = fn X ToFDM 


0.003 and 0.018 


Useful bit rates of DVB-T system 


4.98 Mbits/s for QPSK 
9.95 Mbits/s for 16QAM 


Useful bit rates of 2D LP OFDM 


5.33 Mbits/s for L = 16 


for QPSK 


5.51 Mbits/s for L = 32 
5.60 Mbits/s for L = 64 


Useful bit rates of 2D LP OFDM 


10.66 Mbits/s for L = 16 


for 16QAM 


11.02 Mbits/s for L = 32 
11.20 Mbits/s for L = 64 



Fig. 4 and Fig. 5 give the BER measured at the output 
of the Viterbi decoder for a mobile speed of 20 km/h and 120 
km/h respectively. Note that the value of the pilot power Pp 
has been optimized through simulation search in order to ob- 
tain the lowest BER for a given signal to noise ratio (SNR). 
The performance of the DVB-T system is given with perfect 
channel estimation, taking into account the power loss due to 
the amount of energy spent for the pilot subcarriers. It appears 
in Fig. 4, for low-speed scenario, that the system performance 
is similar to that of the DVB-T system with perfect channel 
estimation. This is due to the power loss due to the pilot which 
is lower with the proposed system. In Fig. 5, for high-speed 
scenario and QPSK, by choosing the spreading lengths offer- 
ing the best performance, there is a loss of less than 1 dB for a 
BER = 10~^, comparing to perfect channel estimation case. 
For 16QAM, the loss is less than 2.5 dB which is really satis- 
fying given that f3 = 0.018, corresponding to a mobile speed 
of 120 km/h. 

5. CONCLUSION 

In this paper, we propose a novel and very simple channel 
estimation for DVB-T. This technique, referred to as spread 
pilot channel estimation, allows to reduce the overhead part 




Eb/No (dB) 

Fig. 3. MSE performance obtained with the analytical ex- 
pression and by simulation ; QPSK ; Speeds: 20 km/h and 
120 km/h ; /3 = 0.003 and 0.018 



dedicated to channel estimation. An analytical expression 
of its MSE, which is a function of the autocorrelation of the 
channel, is given. It allows to highlight and understand that 
the choice of the spreading factors has to be made according 
to the channel characteristics. More generally, this estima- 
tion approach provides a good flexibility since it can be opti- 
mized for different mobility scenarios by choosing adequate 
time and frequency spreading factors. 

This work was supported by the European project CELTIC 
B21C ("Broadcast for the 21st Century"). 

APPENDIX 

In this section, a property from the random matrix and 
free probability theories is defined for the computation of the 
SI variance (6). Furthermore, the computation of the autocor- 
relation function of the channel Rrh is carried out. 

Random matrix and free probability theories property 

Let C be a Haar distributed unitary matrix [5] of size 
[L X L]. C = ( Cp, Cm) can be decomposed into a vector Cp 
of size [L X 1] and a matrix C„ of size [L x (L — 1)]. Given 
these assumptions, it is proven in [7] that: 

CuKCu o^Pu {I - Cpcf) (13) 

where a = 1 is the system load and P„ = 1 is the power of 
the interfering users. 

Autocorrelation function of the channel 

The autocorrelation function of the channel writes: 

Rhh (An, Aq) = E {H,n.s [n, q] i?;„ , [n - An, q - Aq]) 

(14) 
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Fig. 4. Performance comparison between the DVB-T sys- 
tem with perfect channel estimation and the proposed 2D LP 
OFDM ; Speed: 20 km/h ; ^ = 0.003 ; L = 64 ; Pp = 7 
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BQPSK ; Spread Pilot: Lf.16 and Lt.4 ; Pp.7 
#16QAM ; DVB-T system with perfect, chan. est. 
016QAM ; Spread Pilot: Lf.4 and Lt.4 ; Pp.3 
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Fig. 5. Performance comparison between the DVB-T sys- 
tem with perfect channel estimation and the proposed 2D LP 
OFDM ; Speed; 120 km/h ; (3 = 0.018 ; L = 16, 32 and 64 



We can express the frequency channel coefficients H^.s [n, q] 
as a function of the channel impulse response (CIR): 

JVfft-I (-^/+") 

H^^s[n,q]= 7m,9 We^''"^?^'' (15) 

k=Q 

where jm,q [k] is the complex amplitude of the fcth sample of 
the CIR during the gth OFDM symbol of the mth frame, and 
A^FFT is the FFT size. Therefore, by injecting (15) in (14), the 
autocorrelation function of the channel can be rewritten as: 

Rhh (An, Aq) = 

A^FFT — 1 A^FFT — 1 . 

IV^ E E E{7™.,[fc]7;;„-Ajfc']}e-^-^'= 

fc=0 fc'=0 

(16) 

Since different taps of the CIR are uncorrected, it comes: 

RHHiAn,Aq) = 

Wfft-I 

According to Jake's model [8], the correlation of the fcth sam- 
ple of the CIR is: 

E {7m,g [k] l*m.q-Aq W} = PkJf) (27r/l5 AgToFDM) (17) 

where pk is the power of the fcth sample of the CIR, Jo (•) the 
zeroth-order Bessel function of the first kind, Jo the maxi- 
mum Doppler frequency and Tqfdm the total OFDM symbol 
duration. Finally, the autocorrelation function of the channel 
(17) can be expressed as: 

Rhh {An, Aq) = 

E Pfee"'^"^'^ Jo (2^/z3AgroFDM) ^^^^ 
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